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ABSTRACT 



a^ 
o 
\o 

Aims. We address the origin and evolutionary status of hot subdwarf stars by studying the optical spectral properties of 58 subdwarf O (sdO) 
stars. Combining them with the results of our previously studied subdwarf B (sdB) stars, we aim at investigating possible evolutionary links. 
Methods. We analyze high-resolution (R> 18000), high-quality optical spectra of sdO stars obtained with the ESO VLT UVES echelle spec- 
trograph in the course of the ESO Supernova la Progenitor Survey (SPY). Effective temperatures, surface gravities, and photospheric helium 
O ' abundances are determined simultaneously by fitting the profiles of hydrogen and helium lines using dedicated synthetic spectra calculated 
from an extensive grid of NLTE model atmospheres. 



Oh 
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CZ3 Results. We find spectroscopic or photometric evidence for cool companions to eight sdO stars as well as a binary consisting of two sdO stars. 

' A clear correlation between helium abundances and the presence of carbon and/or nitrogen lines is found: below solar helium abundance, no 
^ , sdO star shows C or N lines. In contrast, C and/or N lines are present in the spectra of all sdO stars with supersolar helium abundance. We thus 
? use the solar helium abundance to divide our sample into helium-deficient and helium-enriched sdO stars. While helium-deficient sdO stars are 
r> scattered in a wide range of the Teif-log (g)-diagram, most of the helium-enriched sdO stars cluster in a narrow region at temperatures between 

: 40000 and 50000K and gravities between log g=5.5 and 6.0. 
■ - - ' Conclusions. An evolutionary link between sdB stars and sdO stars appears plausible only for the helium-deficient sdO stars. They probably 
have evolved away from the extreme horizontal branch, i.e. they are the likely successors to sdB stars. In contrast, the atmospheric properties 
of helium-enriched sdO stars cannot be explained with canonical single star evolutionary models. Alternative scenarios for both single star 
(late hot flasher) as well as for binary evolution (white dwarf merger; post-RGB evolution) fail to reproduce the observed properties of helium- 
enriched sdO stars in detail. While we regard the post-RGB scenario as inappropriate, the white dwarf merger and the late hot flasher scenarios 
remain viable to explain the origin of helium-enriched sdO stars. 

Key words, binaries; spectroscopic - stars: abundances - stars; atmospheres - stars: fundamental parameters - stars: horizontal branch - stars: 
subdwarfs 



1. Introduction 

Hot subluminous stars are an important population of faint blue 
stars at high galactic latitudes closely related to the horizon- 
tal branch. A proper spectral classification of hot subluminous 
stars is rendered difficult by the diversity of the helium line 
spectra. They can roughly be grouped into the cooler sdB stars. 

Send offprint requests to: U. Heber 
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* Based on observations collected at the Paranal Observatory of the 
European Southern Observatory for program No. 165.H-0588(A) and 
167.D-0407(A). 



whose spectra typically display no or only weak helium lines, 
and the hotter sdO stars, which have a higher helium abundance 
on average and can even be dominated by helium. The former 
have recently been studied extensively because they are suf- 
ficiently common to account for the UV excess observed in 
early-type galaxie s. Pulsating sdB stars are important tools for 
asteroseismology (Charpinet et al. 2004) and sdB star s in close 
binaiies may qualify as Supernova la progenitors ( Maxted et alJ 
l2000HGeieret al.l2006h . 

Subluminous B stars have been identified as Extreme 
Horizontal Branch (EHB) stars (iHeber-1986.) . i.e. they are core 
helium burning stars with hydrogen envelopes too thin to sus- 
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tain hydrogen burning (unlike normal HB stars). Therefore they 
evolve directly to the white dwarf cooling sequence avoid- 
ing the Asymptotic Giant Branch (AGB). While the sdB stars 
spectroscopically form a homogenous class , a larg e vari ety of 
spectra is observe d among sdO stars llHeber&.Teffervill992t 
iHeher etalJl20n4 . Most subluminous B stars are helium poor, 
whereas only a relatively small fraction of sdO stars are. 

Ev er since the pioneering work by lOreenstein & Sargenll 
lll974 . the helium-rich sdO stars were believed to be linked 
to the evolution of the hydrogen-rich subluminous B stars. Any 
evolutionary link between subluminous B and O stars, how- 
ever, is difficult to explain since the physical processes driv- 
ing a transformation of a hydrogen-rich star into a helium-rich 
one remain obscure. The convectiv e transformatio n has been 
explor ed by IWesemael et alJ (Il982h as well as by Groth et alJ 
(Il985h . While the former found helium convection to occur 
even at subsolar helium abundances which mixes helium from 
deeper layers into the photosphere, the latter concluded that a 
helium driven convection zone develops only in helium-rich 
atmospheres. If the latter is true, convective transformation 
would not work. Non-standard evolutionary models were in- 
troduced to explain the formation of sdO stars (e.g. lSwei^artl 
Il997bt iBrownet alJboOll: iMoehler et al]l2004 . In particular, 
the late hot flasher scenario predicts that the core helium 
flash may occur when the star has already left the red giant 
branch (RGB) and is approaching the white dwarf cooling se- 
quence (delayed He core flash). During the flash. He and C are 
dredged-up to the surface. Hydrogen is mixed into deeper lay- 
ers and burnt. The remnant is found to lie close to the helium 
main sequence, i.e. at the very end of the theoretical extreme 
horizontal branch. 

The fraction of sdB stars in short period b inarie s (peri- 
ods less than ten days) is high. iMaxted et alJ ( I2OOII) found 
2/3 of their sdB sample (drawn from the Palomar Green sur- 
vey) to be such binaries, whereas somewhat lower fractions 
of 4 0% and 48% were found recently for the SPY sample 
jNapiwotzki et al. 200^ and for a sample drawn from the 
Edinburgh-Cape catalog ( Morales-Rueda et al.! '2006). respec- 
tively. Obviously, binary evolution plays an important role for 
the formation of sdB stars and possibly als o for that of the 
sdO stars. A recent population synthesis study jHan et alJ2Q03h 
identified three channels to form sdB stars: (i) one or two 
phases of common envelope evolution, (ii) stable Roche lobe 
overflow, and (iii) the merger of two helium-core white dwarfs. 
The latter could explain the population of single stars. The 
simulations by Han et al. (2003) cover th e observed parameter 
range of sdB stars (see Lisker et al.'2005', henceforth paper I). 

The ESO Supernova la Progenitor SurveY (SPY, 
[Napiwotzki et al. 2001) has obtained VLT/UVES spectra 
for over 1000 white dwarf candidates to test possible scenarios 
for type la supernovae by searching for double degenerate 
white dwarf binary systems close to the Chandrasekhar mass 
limit. Many of the target sta rs of SPY came from the Hamburg 
ESO survey tWisotzki etalJil996) . SPY also observed 137 hot 
subluminous stars that entered the target sample because they 
were previously classified as white dwarfs. 76 of these stars 
were now classified as sdB/sdOB, a nd 58 as O-type subdwarfs 
llChristheb et alJ200ltlLiskeJ2003l) . 



The data for our hot subdwarf sample are of unprecedented 
quality and homogeneity. Spectral analyses of all sdB stars 
from that sample have already been presented in paper I in or- 
der to test evolutionary mo dels, i n particular the binary popu- 
lation models of faan et alJ (120031) . Two diagnostic tools - the 
effective temperature (Teff) vs. surface gravity (log {g)) diagram 
and the cumulative luminosity functi on - yielded conflicting re- 
sults. Moreover, the models of lHan et alJ (f2003h predicted some 
EHB stars to be hotter than the sdB stars contained in the sam- 
ple of paper I. This led us in paper I to the conclusion that their 
sample of hot subdwarfs may not be sufficiently complete to 
describe the whole parameter range covered by the simulations, 
and that it needs to be complemented by a similar analysis of 
subdwarf O stars. 

In this paper we focus on the subdwarf O spectra from the 
SPY sample and present the spectral analysis of the high reso- 
lution spectra using state-of-the-art NLTE model atmospheres. 
This enables us to address the still open question about the evo- 
lutionary status of hot subdwarfs. We combine our results for 
the SPY sdO stars with those for the SPY sdB stars from paper I 
to cover the entire stellar atmospheric parameter range with the 
same high accuracy. This allows for the first time a detailed 
comparison with evolutionary model predictions, particularly 
in the regime of hotter temperatures and higher helium abun- 
dances than in previous studies of sdB stars alone. 

The paper is organised as follows. In section |2] we out- 
line the available data and explain our criteria for spectral 
classification. Section |3l briefly describes the construction of 
a new NLTE model atmosphere grid. In sectionl^we detail the 
derivation of stellar parameters by fitting model atmospheres 
to observed spectral lines and present the results in section [S] 
Various evolutionary scenarios are tested in section |6l leading 
to a discussion and summary in section^ 

2. Observation, data reduction and spectral 
classification 

2.1. Observation and data reduction 

Observations were obtained at the ESO Very Large Telescope 
with UT2 (Kueyen) equipped with the UV- Visual Echelle 
Spectrograph (UVES). A slit width of 2'.'1 was used, re- 
sulting in a resolving power of 18 500 (spectral resolution 
of 0.36 A at Ha) or better Wavelength coverage of 3300- 
6650 A is achieved, with gaps at 4500-4600 A and 5600- 
5700 A JNapiwotzki et alJ200ltfKoester et alJ200lh . For most 
of the stars, two exposures in different nights were taken, since 
SPY was originally intended to search for radial velocity (RV) 
variable objects. The spectra were then reduced with a proce- 
dure developed by C. Karl using the ESO MIDAS software 
package, partly based on the UVES pipeline developed at ESO. 

2.2. Spectral classification 

For many years it has been att empted to establish a consistent 
classificati on of hot s ubdwarfs jGreen et alJ l986 : Jeff erv et alJ 
,1997:.Drilling et alJ2003) that applies to lower resolution spec- 
tra. A transfer of these schemes to spectra of higher resolution 
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would be desirable. However, in case of our UVES spectra, we 
lack the He ii 4542A line due to a wavelength gap between 
the blue and the red arm of the spectrograph. Since this line 
plays a crucial role in the classification scheme of Drillina et al. 
(I2OO3I) . we cannot apply this scheme to our data. Instead we 
aim for a less detailed classification strategy. 

In the process of classifying all SPY spectra llLiskeil2003l) . 
we therefore just separated the spectra into hot white dwarfs 
(of various subtypes), sdB, He-sdB, sdO, and He-sdO stars by 
visual inspection and comparison with spectra of prototypi- 
cal stars. The terms He-sdB/He-sdO were introduced to mark 
extrem ely helium-rich sdB/sdO stars (* Ahmad & .TeffervlbOQSl 
I2OO5I) . Note that we did not use the intermediate class of sdOB 
stars introduced by iBaschek & NorrisI (Il975h - we subsumed 
them in the term sdB. For our programme stars in the present 
paper we continue to use only the general term sdO , in order to 
prevent any bias from qualitative visual classification entering 
our quantitative study of atmospheric parameters. 

It is wor t hwhi le to note that our classification and the one of 
iGreen et'aP ([1986) do not agree as becomes apparent from sev- 
eral stars that we have in common (see Table[2of this paper and 
Table 1 of paper I). Seven programme stars (four s dB and three 
sdO) h ave been misclassified as white dwarfs bv iGreen et alJ 
While none of the sdB stars in paper I has been clas- 
sified as sdO star in the PG survey, six of our sdO stars are 
classified as sdB in the PG catalog. Five of them turn out to 
be helium deficient (see section|5}- It is likely that the helium 
deficient sdO stars are subsumed in the sdB class of the PG 
catalog. Hen ce the fraction of sdO stars is underestimated by 
iGreen et alJ ( fl986i) . 

The high spectral resolution of the UVES spectra does not 
only allow to identify and measure hydrogen and helium lines, 
but also the metal line spectrum can be investigated. Lines from 
highly ionized carbon and nitrogen are found in normal O- 
type stars. For a refined spectral classification we thus used 
the presence or absence of characteristic absorption lines of 
carbon and nitrogen to divide our sample into several sub- 
classes ('CN-scheme'). For carbon we used the strongest lines 
of C HI (4650A, multiplet 1, 4070A, mult. 16, and 4186.9A of 
mult. 18) and of C iv (5801. 33A & 5811.98A, mult. 1, and 
4658. 3A, mult. 8). For nitrogen the strongest lines of N iii are 
used (4640A, mult. 2 and 4858A, mult. 9, 4379. 1 1 A, mult. 18,) 

We adopt the following classification scheme: 

- Stars with C iii and/or C iv visible in their spectra but with 
no N III are classified as 'C. 

- Stars with N iii and C iii and/or C iv in their spectra were 
classified as 'CN'. 

- The absence of C iii and C iv and the presence of N iii leads 
to the classification 'N'. 

- With no lines of C iii, C iv, and N iii visible a star was 
classified as '0'. 

This is not an optimized, quantitative classification scheme, but 
it is intended to demonstrate that even such a qualitative subdi- 
vision already shows a correlation with the atmospheric param- 
eters of our stars, as we shall outline below. The resulting clas- 
sifications are presented in Table ^ We checked whether the 
presence or absence of CN might be significantly influenced 



by the signal-to-noise ratio (S/N), and found that the S/N dis- 
tribution of our spectra does not differ for different CN-classes 
at low and intermediate S/N values, for which it could be cru- 
cial. Hence the CN classification scheme is not significantly 
affected by noise. 

2.3. Search for companions of sdO stars 

In paper I we found a significant fraction of sdB stars with cool 
companions (24 out of 76 sdB stars). Therefore we searched for 
spectroscopic and photometric evidence for such companions 
amongst our sdO stars. 

Only for one subluminous O star, HE 1502-1019 (ahas 
PG 1502-103; EC 15026-1019), spectroscopic signatures for 
a cool companion were found , like the the G-band, and the Mg- 
triplet. iFerguson et alJ ( 1 19841) classified the companion star as 
spectral type K0.5. The weakness of the helium lines implies 
a very low helium abundance of about He/H=0.001. We esti- 
mate the effective temperature from a comparison to synthetic 
spectra (see section|3|i to be near 45 OOOK . 

HE 0301-3039 turns out to be radial velocity variable and 
its spectrum also is composite; it indicates that HE 0301-3039 
is a binary consisting of two sdO stars with spectra dominated 
by helium lines, the first such system found jLisker et all2004l) . 

HE 1200-1924 (ahas EC 12001-1924; Feige 54) shows 
a helium dominated absorption line spectrum typical of a 
sdO star However, Balmer emission lines are superimposed - 
strong and broad Ha, strong but narrow H/3 and weak Hy and 
H6 emission. These emissions do not originate from the photo- 
sphere of tiie sdO star but may stem from a companion star No 
spectral absorptions typical for a cool companion can be found. 
The B - J and J - K colours, however, are unusual for a sdO 
star (see below). 

For the above objects too few spectra are at hand to disen- 
tangle the individual spectra and the quantitative spectral anal- 
ysis of these objects has to be postponed until they become 
available. 

We also searched for photometric evidence for companions 
by complementin g the B-band fluxes w ith IR measurements 
from the 2MASS dSkrutskie et alJ2 006) and DENIS' catalogs. 
Photoelectric B magnitudes (accurate to better than + 0"05) 
were taken from the subdwarf data base ( 0stense nl2OO6l) . If the 
former were unavaila ble, B magnitudes are from the Hamburg- 
ESO survey d Wisotzki et al. 2000. accurate to +0?'2). 

J magnitudes and corresponding errors are provided by 
2MASS for 51 stars, and by DENIS for one star Only the 
brightest stars were detected in the K band. A color-magnitude 
diagram is displayed in Fig.^ Interstellar exti nction and red- 
dening have been corrected for using the maps of lSchlegel et alJ 
til998,) . As most of the stars are located at high galactic lat- 
itudes, ^(B - V) is small (< 0^12) except for the low lati- 
tude star HZ 1. For the latter, the reddening as inferred from 
the Schlegel maps is overestimated, since the star lies near the 
Galactic plane (za!80pc). Instead, we compare the measured 
B -V X.O model predictions, yielding ^(B - V) = 0^25. 



' http;//cdsweb.u-strasbg.fr/denis.html 
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Fig. 1. Colour-magnitude diagram, showing B magnitude ver- 
sus optical-to-near-infrared colour (B-J). Numbers denote ob- 
jects that lie off the bulk of 'normal' sdOs, implying a cool 
companion. The sdO-sdO binary HE030 1-3039 is shown as 
open circle. HE 1502-1019 and HE 1200-1924 are composite 
spectrum objects (see text). 

In paper I we found that sdB stars with spectroscopic sig- 
natures of cool companions are redder than B - J x -0?'4. 
Similarly, eight sdO stars in Fig. 1 lie at B - 7 > -0^5 and 
clearly separate from the bulk of sdOs. We number those ob- 
jects in the figure. Among them is HE 1200-1924, whose emis- 
sion lines provide additional support for the existence of a cool 
companion to that sdO star. 

Only for 12 stars, K band mag nitudes to better than +0^2 
are available. IStark & W ade (2003) studied the distributions of 
J - K colours of sdB and sdO stars and find both to be double- 
peaked, the sdO stars at J - K - -(J?2 and +0^2, respectively, 
indicating that the latter sdO stars have cool companions. J -K 
measurements are available for five of the eight stars high- 
lighted in Fig. n All five have positive J - K indicating the 
presence of cool companions. 

Judged from our experience with sdB stars, we can assume 
the blue spectra of these stars not to be affected by light from 
the cool companion. Those stars are marked in Tableland the 
quantitative spectral analysis is performed in the same way as 
for the other stars (see section|4j. 

3. Stellar atmosphere modeling 

High quality optical spectra like those from SPY demand accu- 
rate theoretical counterparts as input for fitting routines. 

Only few detailed quantitative spectral analyses have been 
carried out for the sdO class because of the complexity of 
their spectra. In addition deviations from local thermodynamic 
equilibrium (NLTE) have to be taken into account, because of 



their high effective temperatures, whereas the LTE assump- 
tion is reasonable for sdB stars. Therefore, o nly few quan- 
titative spectral analyses have been published dPreizler et alJ 
[l990; Thejll et al. 1994), with some conflicting results becom- 
ing apparent. The reliability of the NLTE calculation depends 
strongly on the quality and sophistication of the model atom 
and the atomic data used. The helium model atoms are of ut- 
most importance for the modeling of sdO atmospheres. We 
used more detailed model atoms for He i and He ii than in 
^reizler et al. (1990) and constructed a new grid of atmo- 
spheric models and synthetic spectra using a state-of-the-art 
NLTE model atmosphere code. 

An extensive grid of NLTE atmosphere models was 
calculated using the la test version of the PR02 code 
jWerner & Drei zlei*1999) that employs a new temperature cor- 
rection technique (Dreizler 2003). A new detailed model atom 
for helium appropriate for the sdO temperature regime was 
constructed. 2700 partially line blanketed NLTE model atmo- 
spheres consisting of hydrogen and helium were calculated re- 
sulting in a grid of unprecedented coverage and resolution, 
extending from 30000 K to 100000 K in Teff. The gravity 
ranges from logg=4.8 to 6.4 and the helium abundance from 
log Nhc/Nh^-A- to -1-3 in order to match the diversity of ob- 
served spectra. The step sizes are 2000 K from 30000 K to 
52 000 K and 5 000 K from 55 000 K to 100000 K; 0.2 and 
~0.5 dex, respectively. 

As a test we carried out the spectral analyses of several pro- 
gramme stars using both the old grid of Dreizler et al. ( 1923) 
and the new one. The synthetic spectra from the new models 
match the observed profiles much better in all cases. 

The differences of the current model grid compared to the 
one used by IPreizler et alJ (Il99dl) can be traced back to two 
major improvements: The resonance lines of Hell could only 
be treated in detailed balance in the previous model grid. This 
eliminates the transition from the statistical equilibrium equa- 
tions resulting in a more stable numerical behavior. With an 
improved numerical stability of the code this approximation is 
no longer necessary. Since the resonance lines of He II have 
a strong influence on the structure of the atmosphere of sdO 
stars, this improve ment is the do minant effect. The second 
improvement since IPreizler et alJ (119901) is the treatment o f 
the level dissociation according toMu mmer & Mihalasl l ll988l) . 
This mainly affects higher line series members and is therefore 
important for a precise gravity determination. 



4. Spectral analysis: techniques 

4.1. Line profile fitting 

Atmospheric parameters (Teff, log (g), log (y), y=He/H by num- 
ber) were determined simultaneously by fitting the synthetic 
spectra to obser ved hydrogen and hel ium line spectra using a 
-procedure (Napi wotzki et alJl999l) . 
For each line, the continuum level is determined and nor- 
malized to 1 , in order to compare it to the synthetic spectrum. 
For sdO stars with strong hydrogen Balmer lines, the Ho- was 
excluded in the parameter determination because peculiarities 
in the Hor line core were often reported, perhaps caused by 
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Fig. 2. Sample fit of two sdO star spectra (grey) using the new 
grid of synthetic NLTE spectra as described in Section |3] In 
HE 1135-1134 (top panel) hydrogen is more abundant than 
helium, whereas in HE 0952+0227 (bottom panel) hydrogen is 
a trace element only. The model spectra yielding the best fits 
are shown as black lines. For details see text. 



stellar winds dHeher et alJl2003bh . Nevertheless Ha was kept 
in the final plot for visually examining any significant devi- 
ations from the model profiles. Those sdO stars showing no 
evidence for hydrogen from the visual inspection of the He ii 
Pickering decrement required a modified fitting strategy. If hy- 
drogen were present as a trace element it could be best de- 
tectable as a contribution of Ha to the He ii Pickering line 
6560A but may be too weak to be visible to the eye. In a first 



step Jeff and log (g) are derived by a fit using all lines including 
6560A. Then a second fit is carried out to derive the hydrogen 
abundance keeping Teff and log (g) fixed from the first run. In 
a third step new values for Tetj and log (g) are derived while 
keeping the hydrogen abundance fixed. If the new values differ 
by more than 0.5%, we repeat the procedure until convergence 
is achieved. 

Two stars in our sample ( DeHt 2 and K2-2) have already 
been analysed b v'Napi w otzkil (11999 ) with similar NLTE model 
atmosphere techniques as presented here, and are known to 
be central stars of planetary nebulae (CSPN). DeHt 2 was 
classified as a high luminosity CSPN and Teff =117000 K, 
log(§)=5.64 and log(y)=-0.79 were derived, whereas K 2- 
2 was classified as a peculiar CSPN and reff=67000 K, 
log(^)=6.09 and log (3')=-!. 55 resulted. These stars will not 
be discussed further in this paper 

For 46 out of the remaining 53 stars the fit procedure con- 
verged to convincing line profile fits (as displayed in Fig.|2lfor 
two typical cases). As can be seen the agreement between the 
best fit synthetic spectra and the observed ones is excellent. 

However, for seven stars the fit procedure failed due to var- 
ious reasons. We shall discuss those stars in some detail below. 

HE 1349-2320 (ahas EC 13494-2320) and 
HE 1355-0622 (PG 1355-064): weak He I lines and 
strong He n lines indicate high temperatures. However, Balmer 
and He 11 lines cannot be matched simultaneously. The Ha 
lines display weak emission cores not reproduced by our 
models. 

HE 1512-0331 (see IChristheb et alJl200ll) displays very 
broad Balmer and He 11 lines with emission cores of He 11, 
6560A and 4686A. Its gravity is too high for our model 
grid. Extrapolation suggests that its temperature may be near 
80000 K and log (g) slightly below 7. We regard the star as a 
DAO white dwarf. 

HE 1518-0948 (ahas PG 1518-098, see iMoehler etall 

I1990I) : Weak He I lines indicate high temperature, while 
Balmer and He n lines are narrow but strong indicating that 
it is a hot, helium-rich sdO star of low gravity. In fact, the at- 
tempts to fit its spectrum showed that the gravity must be lower 
than the lower limit of the grid, i.e. log (g)-4.S at a T^ff near 
60000 K. 

HE 2305-1155 displays strong and broad He I lines in- 
dicating a low Teff. Trials to match its spectrum indicated a 
very high gravity beyond the limit of our model grid (i.e. 
log (g) >6.4). 

HE 1008-179 (see IChristheb et alJbOOlT) is a very hot, 
high gravity sdO star showing a Ha emission core. Individual 
Balmer lines can not be matched simultaneously (the so-called 
Balmer line problem, see Napiwotzki 1999). Attempts to fit its 
spectrum indicated that its parameters are probably beyond the 
high temperature, high gravity limit of the model grid, i.e. ex- 
ceeding Teff^lOOOOO K and log (g)=6A, and the star should be 
classified as an extremely hot white dwarf. 

EC 1 1 481-2 303 has already been been analysed by 
IStvs etal.l il200nh who derived reff=41790 K, log(g)=5.84, 
log (37)=-!. 85. The authors were unable to match the UV spec- 
tral energy distribution. The Balmer line problem is obvious in 
our optical spectra from the fit, which indicated a higher tem- 
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Fig. 3. Distribution of relative errors in atmospheric parame- 
ters, as derived from two or more exposures. From top to bot- 
tom the distributions for T^ff, g and y - Nhc/Nh are shown. 



and the low noise level of our data. Nevertheless we decided 
to derive a more reliable estimate of the true errors on the ba- 
sis of the existing two or more exposures of each star, where 
available. We determined the relative error for T^ff, g and y 
by taking the difference between the fit results of each expo- 
sure and dividing by the weighted mean of these fit results (see 
Section im . Figure |3l shows the histogram distribution of the 
resulting relative errors. We find the l-cr- value (sample stan- 
dard deviation) of these distributions to yield global relative 
eiTors of Areff/reif = 0.026, Ag/g = 0.25 and Ay/y = 0.30. 
This corresponds to AlogTeff = 0.011, Alog(g) = 0.097 and 
Alog(3;) - 0.11. We note, that some relative errors are par- 
ticularly large, i.e. more than 2-cr (sample standard deviation) 
of these distributions. These so-called "outliers" are flagged in 
Tab.n 



5. Spectral analysis: Results, trends and biases 



As outlined in lEdelmann et alJ (l2003l) and paper I we shall 
search for trends in parameter space. We shall discuss the distri- 
butions of gravity, temperature, and helium abundance, taking 
into account the CN-classification scheme (see section lT^ . As 
systematic differences in atmospheric parameters lead to selec- 
tion biases, we discuss this issue at the end of this section. 



perature 50000 K), low er helium content («; -2.8) but sim- 
ilar gravity than found bv lStvs et al.l(l2000l) . The Balmer line 
problem in ve ry hot wh ite dwarfs was traced back to metal line 
blanketing by IWerneJ (il996) . .Hammer et al.l (l2003l) analysed 
the UV spectrum of EC 11481-2303 and derived very large 
iron and nickel abundances of 10 to 100 times solar, causing a 
strong line blanketing effect. Hence, the neglect of metal line 
blanketing is the likely cause of our failure to match the optical 
spectra of this star 

Because of the peculiarities described and problems en- 
countered in the fit procedure, we do not include these stars 
in the further discussion and are left with a working sample of 
46 stars. 

Since the SPY project was originally intended to search for 
radial velocity variable stars, two or more exposures of each 
star, if available. Thus the final values were calculated as the 
average of the individual fit results, weighted with the S /N of 
the corresponding spectra. The final values are presented in 
Table [2 We note, that for seven stars of the working sample 
only a single (useful) exposure was available and denote them 
by 'S' in Table [2 Effective temperatures range from 36 000 K 
to 78 000 K, gravities from log(g)-4.9 to 6.4, and helium-to- 
hydrogen ratios from 0.0003 to more than 1000. 

By adopting a canonical mass of 0.47 Mq for our pro- 
gramme stars, we can further derive the luminosity in solar 
units (Table 0. 

4.2. Error determination 

The statistical 1-cr-errors from the fit procedure are typically 
lower than 100 K, 0.04 dex and 0.04 dex for T^s, \ogig) and 
logiy), respectively, as a consequence of the high resolution 



5.1. Atmospheric parameters and C & N line strengtlis 

The following conclusions can be drawn from the distribution 
of stars in the Teff-log (y)-diagram (Fig.|3: 

(i) SdO stars with a hehum abundance below solar, ranging 
from 0.0003 < He/H < 0.08, are scattered in a wide range 
of effective temperatures from 35 kK ^ T^ff $ 80 kK, and no 
correlation can be found. 

(ii) SdO stars with a helium abundance exceeding the solar 
value, ranging from 0.1 5 He/H S 1, tightly cluster around 
Teff ~ 40 kK and He/H!^0.5. 

(iii) Above a helium abundance of He/H=10 the sdO stars 
show a clear dependence of the helium abundance on tempera- 
ture: the helium abundance decreases with increasing Teff , op- 
posite to the general trend for sdB stars. 

Most strikingly a clear correlation between helium abun- 
dance and CN class becomes apparent from Fig.|3 None of the 
sdO stars with subsolar helium content shows carbon and/or 
nitrogen lines and therefore all of them are of CN class 0. The 
opposite is true for sdO stars with supersolar heUum content - 
all of them show carbon and/or nitrogen lines and therefore are 
either of class C, N, or CN. 

This suggests that the sdO stars should be grouped accord- 
ing to helium content into two classes. Those with supersolar 
helium abundances will further on be referred to as helium- 
enriched sdO stars while those with subsolar helium abun- 
dances will be termed helium-deficient sdO stars. 

Figs. |5] and |6lshow the distribution of the full SPY sam- 
ple of sdO stars in the Teff-log (g)-plane, in comparison with 
the sdB stars f rom paper I. The location of the EHB band 
dPorman et alJl993 ) is shown, along with the helium zero-age 
main sequence (He-ZAMS) which indicates the boundary be- 
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Table 1. Results of the quantitative spectral analysis of the ESO SPY sdO stars sorted by increasing helium abundance (log(y)). 
The results are mean values from two or more exposures, if available. "Outliers", i.e. values with relative errors between two or 
more exposures above the 2-cr level of the global relative error, are flagged by an asterisk (see section]^}. The CN-classification 
is given as des cribed in Section 12. 2 | B magnitudes are from the Hamburg ESO objective prism survey (photographic, accurate to 



+0.2 mag, see lwisotzki if not specified otherwise. 



Object 


ICRS 2000.0 


B 






log (g) 


log (y) 


log(L) 


CN 


Rem. 






mag 


mag 




r -1-1 

[cm s "J 




FT 1 






helium deficient: 




















HE 1059-2735 


11:01:24.95 -27:51:42.9 


15.21" 


3.11 


40966 


5.38 


-3.62 


2.16 





2 


HE 1130-0620 


11:32:41.60 -06:36:54.4 


15.76" 


4.05 


48122 


5.84 


-3.07 


1.98 





1 


HE 1237- 1408 


12:39:56.64 -14:24:48.4 


15.97" 


3.02 


39683 


5.33 


-2.99 


2.15 





1,2,5,8 


HE I3I8-21I1 


13:21:15.66 -21:27:18.5 


14.48" 


3.31 


36254* 


5.42 


-2.91 


1.92 





2,5 


HE 1 115-0631 


11:18:11.69 -06:47:33.2 


14.77" 


4.17 


40443 


5.80 


-2.59 


1.72 





1 


PG 0026+136 


00:28:52.33 +13:54:45.8 


15.74' 


3.17 


38830 


5.38 


-2.39 


2.07 





1 


EC 09445-0905 


09:47:03.39 -09:19:50.5 


15.69" 


4.59 


73862 


6.22 


-2.08 


2.28 





2,4 


HE 1423-0119 


14:25:51.29 -01:33:17.4 


16.56" 


3.12 


52662 


5.50 


-1.61 


2.48 





1 


HE 15 1 3-0432 


15:16:19.17 -04:43:58.0 


16.04" 


4.41 


42699 


5.92 


-1.36* 


1.69 





1,7 


HE 0040-4838 


00:42:31.08 -48:22:16.2 


16.07 


3.79 


41823 


5.67 


-1.35 


1.91 







HE 1047-0637 


10:50:28.79 -06:53:25.9 


14.42" 


1.49 


60650 


5.03 


-1.34 


3.39 





1,5 


HE I356-I6I3 


13:59:12.52 -16:28:01.8 


16.18 


4.11 


55925* 


5.93 


-1.30 


2.15 





5 


HE 1106-0942 


11:09:08.22 -09:58:48.6 


16.34 


4.92 


79742* 


6.40 


-1.03 


2.27 





1 


helium enriched: 




















HE 1238-1745 


12:41:01.16 -18:01:59.0 


14.31" 


3.83 


38219 


5.64 


-0.55 


1.78 


N 


2,5 


HE 1258+0113 


13:00:59.23 +00:57:11.7 


16.23" 


3.80 


39359 


5.64 


-0.53 


1.83 


N 


1,5,7 


HE 2218-2026 


22:21:13.02 -20:11:17.6 


16.28 


4.40 


38330 


5.87 


-0.35 


1.56 


CN 




HE 1135-1134 


11:38:10.66 -11:51:03.8 


15.45 


3.88 


40079 


5.68 


-0.26 


1.82 


N 


1 


HE 1136-2504 


11:39:10.21 -25:20:55.5 


13.83" 


4.25 


41381 


5.84 


-0.16 


1.72 


N 


2,5 


HE 1310-2733 


13:12:50.65 -27:49:03.1 


14.38 


3.76 


40000 


5.63 


-0.15 


1.87 


N 


2,5 


HE 2359-2844 


00:01:38.48 -28:27:42.8 


15.7 


3.85 


38325 


5.65 


-0.15 


1.77 


CN 


3 


PG 2204+070 


22:07:16.20 +07:18:36.0 


15.74" 


3.72 


40606 


5.62 


0.07 


1.90 


N 


1,4,8 


HE 1256-2738 


12:59:01.48 -27:54:19.3 


16.29" 


4.04 


40029 


5.68 


0.09 


1.82 


CN 


2,5 


HE 2203-2210 


22:06:29.38 -21:56:00.0 


15.04 


4.95 


47343 


6.20 


0.45 


1.59 


CN 




HE 1142-2311 


11:44:50.15 -23:28:18.0 


15.37" 


3.86 


54718 


5.80 


0.68 


2.21 


c 


2,5,6 


HEOlll-1526 


01:13:38.24 -15:11:02.6 


14.59" 


3.86 


39152 


6.31 


0.83 


1.16 


CN 


3,8 


HE 1251+0159 


12:54:08.35 +01:43:24.0 


15.24" 


4.47 


48208 


5.98 


1.03 


1.80 


C 


1 


HE 1511-1103 


15:14:17.04 -11:14:13.6 


14.78" 


3.88 


42298 


5.68 


1.10 


1.91 


CN 


1,5 


HE 1430-0815 


14:33:36.93 -08:28:24.8 


15.69" 


2.32 


61011 


5.26 


1.17 


2.97 


C 


1,2,7,8 


HE 1203-1048 


12:05:56.59 -11:05:29.4 


15.69 


3.98 


45439 


5.91* 


1.36* 


1.91 


c 


1 


HE 1446-1058 


14:49:24.49 -11:11:19.0 


15.00" 


3.91 


45000 


5.76 


1.37 


1.94 


CN 


2 


HE 0342-1702 


03:44:58.82 -16:52:42.2 


14.75 


4.18 


41914 


5.78 


1.40 


1.79 


N 




HE 0952+0227 


09:55:34.57 +02:12:47.9 


14.72 


3.84 


46524* 


5.75 


1.41 


2.00 


C 


1,5 


HE 2347-4130 


23:50:19.70 -41:14:01.1 


15.16 


4.07 


44875 


5.83 


1.44 


1.88 


c 




HE 0414-5429 


04:15:30.23 -54:21:58.7 


14.60 


3.86 


44678 


5.76 


1.57* 


1.96 


C 




HE 0914-0314 


09:17:15.62 -03:53:57.3 


14.93 


3.97 


45136 


5.79 


1.65 


1.94 


C 


1 


BPS CS 22955-0024 


20:23:50.26 -25:08:28.8 


15.50'^ 


4.04 


44622 


5.80 


1.72 


1.89 


CN 


4 


HE 0155-3710 


01:58:01.44 -36:56:21.9 


15.08 


4.24 


41405 


5.77 


1.76 


1.76 


N 


3,7 


HE 1136-1641 


11:38:54.62 -16:58:13.4 


14.82" 


4.03 


44646 


5.80 


1.81 


1.88 


C 


2 


HE 0958-1 151 


10:00:42.64 -12:06:00.0 


13.81" 


3.29 


44125 


5.51* 


1.85 


2.18 




1,2,4,5 


HE 0016-3212 


00:18:53.22 -31:56:01.7 


14.27 


3.73 


41674 


5.70* 


2.14 


1.97 


CN 


3,6 


HE 0031-5607 


00:34:07.75 -55:51:05.9 


15.45 


5.37 


41423 


6.25 


2.26* 


1.31 


N 


8 


HE 1316-1834 


13:19:16.93 -18:49:52.0 


16.24 


3.55 


42811 


5.56 


2.32 


2.05 


N 


2,8 


PG 2258+155 


23:00:57.75 +15:48:39.8 


15. 14*^ 


4.77 


42084 


6.08* 


2.84* 


1.56 


N 


1,4 


PG 1632+222 


16:34:16.09 +22:11:40.9 


15.13" 


4.77 


39384 


6.16 


2.84 


1.31 


C 


1,4,8 


HE 0001-2443 


00:04:31.01 -24:26:21.1 


13.66 


4.62 


40975 


5.94 


2.97 


1.60 


N 


3 


HZ 1 


04:50:13.52 +17:42:06.2 


12.60" 


3.95 


41344 


5.68 


3.00 


1.87 


N 


3,4 



1: Also in PG survey (Green et al. 1986) 
3: Also in MCT survey (Demers et al. 1987) 
5: Christlieb et al. (2001) 

7: Possibly cool companion indicated by i? - 7 colour 
° From subdwarf data base (0stensei^OO&) 



2: Also in EC survey (Stobie et al. 1987) 
4:^cCoo^^^ion ^99^, misclassified as white dwarf 
6: Possibly cool companion indicated hy B - J , J - K colours 
8: Only one spectrum available 

* from SDSS g' and r' using the calibration of^mitl^^l^ (200Z). 
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Fig. 4. Helium abundance versus effective temperature for sdO 
as shown by their CN-classification and co mparison to the sdB 
stars (open squares) from the SPY project jLisker et alJ2 005). 
Open circles denote stars of CN class 0, filled triangles are ei- 
ther class N (upside down triangle) or class C, whereas stars of 
CN class are marked by a filled diamond. sdB stars are marked 
by open squares. Measurement uncertainties are given in the 
upper left corner 



5.5 



O 



6 - 



6.5 



I 



O 

o o 



y 

/u 

□ / □ 

□ □ 

a /a ° / 

7 u \' 



O 0.8>~..p / 

o>.^ 



80 70 60 50 40 30 

Tef, [kK] 

Fig. 5. Helium-deficient sdO stars: Distribution of effective 
temperature and surface gravity. The location of the EHB band 
and the helium zero- age main sequence (He-ZAMS, labeled 
by stellar mass in solar units) are also indicated. Notation as in 
Fig. 0] Measurement uncertainties are shown in the upper left 
corner 
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Fig. 6. Helium-enriched sdO stars: Distribution of effective 
temperature and surface gravity. The location of the EHB band 
and the helium zero-age main sequence (He-ZAMS, labeled 
by stellar mass in solar units) are also indicated. Notation as in 
Fig.0] Measurement uncertainties are shown in the lower right 
corner. 



yond whic h no s table helium core burning can be established 
i|Paczvnski il97ll) . 

Both helium-deficient (Fig.|5} and helium- enriched (Fig.|6} 
sdO stars are found at higher temperatures than the sdB stars, 
most of them between 37 kK and 47 kK. They clearly lie out- 
side the EHB band, which is defined as the region between the 
zero-age EHB (ZAEHB) and the terminal age E HB (TAEHB ) 
given by the evolutionary calculations of Dorman et alJ 
A significant fraction (7 out of 33) of helium-enriched sdOs 
appear to lie below the helium main sequence. In view of the 
non-Gaussian distribution of errors (see Sect.|5}, their gravities 
might have been overestimated. Indeed, two stars are amongst 
the "outliers" and the results of another three are based on sin- 
gle spectra. It may therefore be premature to draw reliable con- 
clusions. Additional observations are necessary to derive more 
precise gravities. 

While the helium-deficient sdO stars are scattered in a wide 
Teff-log (g)-range, most helium- enriched sdOs populate a rela- 
tively narrow region (T^f^ from ~40 to ~46kK and log (g) from 
-5.5 to -5.9). 

5.2. Selection bias 

The target objects of SPY were selected from a magnitude lim- 
ited list of candidate white dwarfs. Therefore all subdwarfs 
present in the SPY dataset have only been included "by acci- 
dent", because they were classified erroneously as white dwarfs 
mostly from low resolution objective prism spectra. For the sdB 
stars, we found in paper I that any selection effects cannot be 
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significantly different from typical selection effects present in 
earlier studies of sdB stars. Besides the usual biases of magni- 
tude limited surveys there may be another selection effect in- 
troduced by the SPY observing strategy. As the SPY project 
aimed at including a significant fraction of helium-rich white 
dwarfs, preference was given to DB candidates, many of which 
turned out to be sdO stars. There may thus be a bias in favour 
of selecting sdO stars relative to sdB stars in the SPY project. 
This must be taken into account when the combined sdB and 
sdO sample from SPY is compared to predictions from evolu- 
tionary calculations. 



6. Evolutionary status 

As outlined above, the distribution of the programme stars in 
the Teff-log (g)-diagram (see Figs. |5] and |6jl as well as with 
respect to CN class suggests that helium-enriched sdO stars 
form a population different from the helium-deficient sdO stars. 
Published evolutionary scenarios try to explain the origin of 
sdB and sdO stars either by canonical or non-canonical evo- 
lution of single stars or by close binary evolution with mass 
exchange and common envelope episodes. We now test these 
scenarios by comparing our observational results to their pre- 
dictions. 



6.1. Canonical evolution of single stars 



Canonical EHB models (e.g. Dormanet al. 1993) are charac- 
terized by a core mass of nearly half a solar mass and a tiny 
inert hydrogen-rich envelope. The core mass is fixed by the on- 
set of the core helium flash at the tip of the red giant branch 
and depends only slightly on metalicity and helium abundance. 
Hence the canonical core mass is restricted to a very narrow 
range of 0.46 to 0.5 Mq. This configuration prevents an EHB 
star from ascending the AGB. The post-EHB evolution pro- 
ceeds towards higher temperatures until the white dwarf cool- 
ing track is reached and gravity increases. 

The problem for the formation of EHB stars in this case is, 
how almost the entire envelope of the RGB progenitor is lost at 
precisely the same time as the core reaches the mass required 
for the helium flash. Enhanced mass- loss during or after the 
red giant branch has been postulated jD'Cruz et alJI 19961) but 
no physical mechanism has yet been identified. 

In Fig. we compare the distribution of sdO stars to the 
position of the EHB-band and to post-EHB evolutionary tracks 
and find that none of the programme stars lies on the EHB- 
band. However, it is premature to conclude that they have 
evolved from the EHB. Most of the helium-enriched sdO stars 
cluster in a narrow region of the Teff-log (g)-diagram. The cal- 
culations of lDorman et alJ ( fl993 ) indicate that the pace of evo- 
lution does not change very much through post-EHB evolution. 
Hence post-EHB evolution can not explain the clustering of 
helium-enriched sdO stars. Moreover, the post-EHB scenario 
fails to explain the surface enrichment of helium and it also can 
not explain why the C and/or N lines in helium- enriched sdO 
stars are stronger than helium-deficient ones, since no dredge- 
up process is predicted to occur. 
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Fig. 7. Tracks for post-E HB evolution (black dotted lines) 
from lDorman et al1(ll993l) are superposed onto the Tgff-log (^)- 
diagram for helium-deficient and helium-enriched sdO stars, 
with core masses as indicated next to the tracks. Notation as 
in Figs. m and |6l We also show post-AGB tracks (dashed lines) 
from Schonbernen lll979lll983i) with white dwarf (final stage) 
masses as indicated next to the tracks. The helium main se- 
quence is labeled by stellar mass in solar units. 



The helium- deficient sdO stars, however, can well be ex- 
plained by the post-EHB hypothesis as they have similarly low 
helium abundances as sdB stars. In addition they are scattered 
in the Tetj-log (g)-diagram, hence no slow-down in their evolu- 
tion has to be invoked. Therefore it is reasonable that helium- 
deficient sdO stars are post-EHB stars and have evolved from 
sdB stars, while helium-enriched sdO stars have not. 

Some of the more luminous sdO stars may not be related to 
the EHB at all but may have evolved off the AGB. Therefore 
we included post-AGB evolutionary tracks from Schonbe rned 
(Il979l Il983l) in Fig. H The position of the helium-deficient 
sdO HE 1047-0637 (T^eif ~ 60 OOO/T, log {g) « 5.0) is matched 
by these tracks, indicating a possible post-AGB evolutionary 
stage while the helium- enriched sdO HE 1430-0815 (Teff ~ 
61 QQQK,\og{g) ^ 5.3) may be either a post-AGB or a post- 
EHB star. The short evolutionary timescales ( about 30000 
years from the AGB towards a pre-white dwarf, ISchonberneJ 
drastically reduce the probability of finding true 
post-AGB stars. 

6.2. Non-canonical evolution of single stars 

Since canonical single star evolution can not explain the 
helium-enriched sdO stars, in particular their clustering in the 
Teff-log (g)-plane, the helium enrichment and the C and/or N 
fine strengths, other scenarios have to be investigated, such as 
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Fig. 8. Helium-enriched sdO stars: comparison with an evo- 
lutionary track for an EHB star formed by a delayed helium 
flash in the effective temperature vs. surface gravity plane. The 
track settles onto the helium main sequence (labeled by stellar 
mass in solar units). Symbols as in Fig.|6] 



the late hot flasher scenario, in which the core helium flash oc- 
curs when the star has already left the RGB and is approach- 
ing the white dwarf cooling sequence (delayed helium core 
flash). During the flash. He and C is dredged-up to the surface 
(ISweigarJl997ah . 

Fig. |8l compares an evolutionary track for the late hot 
flasher scenario with the distribution of our /zelium-enriched 
sdO stars. The star encountered a late flash on its way from 
the RGB towards the white dwarf regime. This late hot flash 
forces the star to land on or near the helium main sequence, 
i.e. at the extreme end of the EHB. The final composition 
of the envelope is helium dominated with Y=0.814, comple- 
mented by hydrogen (X=0.154) and enriched with carbon at 
Z=0.032 (or nitrogen if the hydrog en burning durin g the he- 
lium flash phase burns '^C into '^N; ISweigartll997al) . Indeed, 
most of our observed helium-enriched sdO stars lie near the 
model track, suggesti ng that they indeed may originate from 
this scenario (see also iLemke et al.lf 19971) . However, the evo- 
lutionary time scales (1.95x10^ yrs for the evolution shown in 
Fig- EI Sweigart 1997a) are much shorter than for the core he- 
lium burning phase. Accordingly the stars should accumulate 
near the end of the track, i.e. near the helium main sequence, 
which is not the case for our programme stars. 

Although the late hot flasher scenario can explain the he- 
lium enrichment and the line strengths of C and/or N lines as 
due to dredge up, it fails to reproduce the distribution of the 
stars in the Teff-log (g)-diagram in detail. 



6.3. Binary evolution 

As the evolutionary scenarios for single stars discussed above 
partially fail to explain the observed properties of helium- 
enriched sdO stars, we now focus on binary evolution. Two 
flavours of close binary evolution have been envisaged to ex- 
plain sdB and sdO stars. The formation of core helium burning 
EHB stars through close binary evolution has extensively been 
investigated by Han et al. ( 2002. 2003) using the binary popu- 
lation synthesis approach. 

However, it may not be taken for granted that hot sublumi- 
nous stars are indeed core helium burning. They may form from 
red giants that left the RGB before igniting helium in the core 
and evolve through the EHB region as helium stars towards the 
white dwarf cooling sequence, see Heber ( 2003). Such progeni- 
tors of heliu r n core w hite dwarfs, indeed, have been discovered. 
iHeber et alJ ll2003al) found that the sdB star HD 188112 has a 
mass of O.23M0, too small to sustain helium burning. 

6.3.1 . Binary evolution scenarios and the role of white 
dwarf mergers 

iHan et al] 12003. hereafter HPMM) showed in their binary pop- 
ulation synthesis study that three channels are relevant for the 
formation of hot subluminous stars involving either common- 
envelope ejection, stable Roche lobe overflow or a merger of 
two helium white dwarfs. Stable Roche lobe overflow is pre- 
dicted to lead to composite spectrum systems consisting of a 
hot subdwarf and a main sequence star There are eight such 
binaries (sdO -H main sequence star) present amongst our pro- 
gramme stars. 

In paper I, we compared the atmospheric parameters of the 
SPY sdB stars to the HPMM models by using two diagnostic 
tools, namely the Teff-log (g)-diagram and the cumulative lumi- 
nosity function. Our analysis of sdO stars extends the subdwarf 
sample to higher temperatures which may allow to study the 
link between sdB and sdO stars. One would thus expect a com- 
parison of our full sdB/sdO sample with the HPMM models to 
yield more robust results than in paper I. However, we must 
observe selection biases as discussed in section I5T2I 

We present this comparison for the Teif-log (^)-diagram in 
Fig.m The simulation set No. 10 of HPMM was chosen be- 
cause it came closest to the SPY-sdB distribution (see paper I). 
The grey shading of the rectangular areas corresponds to the 
respective number of simulated stars they contain (cf. paper I). 
Higher number densities of simulated subdwarfs correspond to 
darker grey shading. We refer the reader to paper I for more 
details about the general specifications of the HPMM simula- 
tions. 

From the direct comparison of our derived T^a-logig)- 
values to the HPMM simulations, we see two effects. First, 
sdO stars significantly exceed even the hottest temperatures 
that result from any HPMM simulation set, with stars reaching 
up to 80 kK. Even for the sdO stars with lower temperatures, 
no set is able to reproduce their rather wide range in surface 
gravity. HPMM therefore covers only < 38% of all sdO stars. 
Second, by restricting our analysis to stars which come close 
to the HPMM predictions, i.e. those that are apparently con- 
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nected with the sdB sample, a strong disagreement of the ob- 
servational data with the simulation set becomes obvious: the 
relative amount of hot (sdO) and cool (sdB) stars differs signif- 
icantly. Apart from possible limitations of the HPMM models 
(as discussed below in section^, there may be an observational 
bias in favour of sdO stars, rendering a quantitative compari- 
son with HPMM, as performed for the sdB stars in paper 1, very 
difficult or even impossible. 

Let us now focus on the white dwarf merger channel, which 
is supposed to result in stars with very low hydrogen enve- 
lope mass. Hence stars resulting from mergers of two helium- 
core white dwarfs are expected to be found near the hot end 
of the HPMM distribution. Mergers can also produce stars in 
a wider range of masses than the other channels predict. The 
merging process probably induces lots of mixing of nuclear 
processed material to the stellar surface of the remnant, poten- 
tially leading to a helium- and nitrogen-rich surface composi- 
tion. Depending on the efficiency of nuclear burning and mix- 
ing, carbon may also be enriched at the surface. As these pre- 
dictions match the observed properties of the helium-enriched 
sdO stars at least qualitatively, we regard the merger of helium 
core white dwarfs as a viable scenario. There are two predic- 
tions we can test: (i) Stars formed from a merger should not be 
radial velocity variable (unless they stem from triple systems) 
and (ii) the helium enrichment is accompanied by enrichment 
of nitrogen and/or carbon. While our CN-classification scheme 
(section l2T2t is based on the presence of C and/or N-lines in 
the spectrum, the latter issue requires a quantitative abundance 
analysis. We shall address these tests in forthcoming papers. 
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Fig. 9. Comparison of the atmospheric parameters of helium- 
deficient sdO, helium- enriched sdO, and sdB from t he SPY 
project to simulation set No. 10 of ( Han et al]l2003l) . which 
was found to match the sdB distribution of Paper I best. The 
theoretical predictions are shown as shaded Teff-log (g)-boxes, 
where a higher subdwarf density per box corresponds to darker 
shading. The grey scale is shown below the figures. Notation 
as in Fig.0] 



6.3.2. Non-core helium burning stars in close binary 
systems 

The evolution of RGB stars whose envelopes get almost com- 
pletely stripped by Roche lobe overflow in a close binary sys- 
tem before helium burning starts has been investigated e.g. by 
lDriebeetalJlll998f) . These tracks were calculated from a IM© 
model sequence starting from the pre-main sequence stage up 
through the RGB. Large mass loss rates were then adopted and 
the evolution of the resulting helium star was followed. The 
remnant finally evolves into a helium core white dwarf. 

In Fig. ^1 we compare the pro gramme stars to the predic- 
tions of the evolutionary models o flPriebe et alJ (1 1998*) for dif- 
ferent masses. While the position in the Teff-log (g)-diagram of 
any of our programme stars can be matched by a post-RGB 
track of appropriate mass, it is striking that most of the helium- 
enriched sdO stars agree reasonably well with the theoreti- 
cal predictions for a rather narrow mass range, i.e. between 
0.3Mq and 0.33Mq. Hence the distribution of helium- enriched 
sdO stars could be a sequence of low mass stars evolving into 
helium white dwarfs. It is also surprising that the predicted 
masses are close to the minimum mass for the helium main 
sequence, which might be purely coincidental. 

However, the models of IPriebe et alJ lll998h predict rather 
thick hydrogen layers. Therefore, this scenario has difficulties 
to explain the high helium abundances at the stellar surface as 
observed in helium-enriched sdO stars. During the post-RGB 



evolution hydrogen-shell flashes occur which possib lv may 
lead to a dredge up of helium and nitrogen jDriebe et al..l999.) . 
This may explain nitrogen strong-lined objects. It is, how- 
ever, not evident how ca rbon-strong o bjects could be formed. 
Moreover, the models of Driebe et aP ( il999.) indicate that the 
occurrence of hydrogen flashes is constrained to a lower mass 
range of 0.21 to 0.3Mo. 

If the hydrogen were removed from the envelope by some 
unknown process, the evolution would be speeded up consider- 
ably, because most of the luminosity is provided by the hydro- 
gen burning shell, which would be extinguished if the hydro- 
gen mass is too low. Rapid evolution reduces the detectability 
of such stars. 

Therefore we regard the post-RGB scenario as unlikely. 
However it is a viable scenario to explain some of the helium- 
deficient sdO stars since helium core white dwarf progenitors 
are known albeit ra re amongst the sdB stars (e.g. HD 188112, 
iHeheret alJ2003l. 

7. Summary and Discussion 

We analyzed high resolution optical spectra of 58 sublumi- 
nous O stars. We found spectroscopic and/or photometric evi- 
dence for cool companions to eight subluminous O stars. While 
HE 1502-1019 and HE 1513-0432 display only weak helium 
lines, the others have strong helium Unes. HE 0301-3039 is an 
unique binary consisting of two sdO stars with spectra dom- 
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Fig. 10. Tracks for post-RGB evolution (solid lines) towards 
helium core white dwarfs of 0.3Mq to O.414M0 (Driebe et al. 
Il998l) are superposed onto the Teff-log (g)-diagram of helium- 
enriched sdO stars, with masses of the model star as indicated 
next to the tracks. The helium main sequence is labeled by stel- 
lar mass in solar units. Notation as in Fig.|3] 

inated by helium lines (iLisker et alJl2004l) . Evidence for cool 
companions to the other helium-enriched stars rests on pho- 
tometry only. Near-infrared spectroscopy is required for classi- 
fying them. As the photometric measurements for our sample 
are only available for 52 out of 58 stars, the fraction of sdO 
stars with cool companions may be slightly larger. 

A grid of synthetic H-/He-line spectra calculated from 
NLTE model atmospheres was used to derive the stars' at- 
mospheric parameters. Reliable atmospheric parameters were 
obtained for 13 helium- deficient sdO and 33 helium-enriched 
sdO stars. A clear-cut correlation between CN class and helium 
abundance was found. Helium-deficient sdO stars did not show 
any C and/or N lines, while all helium-enriched sdO stars do 
show either C or N lines or both, suggesting that they form two 
different populations of stars. The helium-deficient sdO stars 
are scattered in the Teff-log (g) diagram, whereas most helium- 
enriched sdO stars cluster in a narrow range. 

Comparing the observed distribution to the predictions of 
evolutionary calculations for single as well as for close binary 
stars, we conclude that many helium- deficient sdO stars can be 
explained as evolved sdB stars. Both classes of star are helium 
deficient and have weak metal lines which is caused by atmo- 
spheric diffusion processes. 

Most of the helium- enriched sdO stars cluster in a nar- 
row region of the Tefj-log (§)-diagram at temperatures between 
40kK and 50kK. While diffusion is probably causing helium 
deficiency, it is unlikely to account for the helium enrichment. 
Non-standard evolutionary scenarios were therefore considered 
as well. The predictions from the late hot flasher scenario as 



well as the helium white dwarf merger scenario are roughly 
consistent with the observed distribution of helium-enriched 
sdO stars but do not match them in detail. The occurrence of 
a delayed helium core flash as well as the merger of two he- 
lium white dwarfs may explain the helium enrichment. In both 
cases carbon and/or nitrogen can be dredged up to the stellar 
surface, which would explain the strength of the C and/or N 
lines in helium- enriched sdO stars. 

Some high gravity helium-enriched sdO stars may lie be- 
low the helium main sequence, which is at variance with any 
core helium burning model. Therefore, we considered mod- 
els for post-RGB stars with inert helium cores which evolve 
through the sdB/sdO regime in the Teff-log (g) diagram into he- 
lium core white dwarfs. This is the only scenario that can ex- 
plain stars to lie below the helium main sequence. However, 
it is not obvious how the helium enrichment is brought about 
as the hydrogen envelopes of post-RGB models has to be rel- 
atively thick or, otherwise, the stars would evolve too fast to 
be observable in large quantities. If the existence of a popula- 
tion of high gravity sdO stars below the helium main sequence 
could be confirmed by surveys with larger sample size, a pos- 
sible explanation could be provided by the post-RBG scenario. 

Our conclusions can be tested by measuring the binary fre- 
quency and the C and N abundances. If helium-deficient sdO 
stars are evolved sdB stars, their binary fractions should be 
the same. As about 40% of the sdB stars are in close b ina- 
ries with periods below 10 days j Napiwotzki et alJl2004h . we 
expect 5 helium-deficient sdO stars in our sample to be radial 
velocity variable on time scales of 10 days or less. An inves- 
tigation of the binary frequency is underway in our group. If 
helium-enriched sdO stars result from mergers, no radial ve- 
locity variations would be expected (except for objects arising 
from triple stars). If material processed by nuclear burning is 
dredged up, the helium enrichment is expected to be accom- 
panied by enrichment of nitrogen and/or carbon. While our 
CN-classification scheme (section lT^ is based on the presence 
of C and/or N-lines in the spectrum, the latter issue requires a 
quantitative abundance analysis to test the late hot flasher and 
the merger scenario. 

In addition, improvements to our analyses of the SPY- 
sample are required for a more detailed comparison with pre- 
dictions of evolutionary calculations. The role of metal-line 
blanketing in the NLTE model atmospheres needs to be in- 
vestigated. Due to the complexity of the problem, only few 
studies of N LTE metal lin e blanke ting in sdO st ar atmospheres 
are available faaas et alJ l996: La nz et alJl997l lbeetien 20o3) 
which indicate that the effective temperatures may have to be 
reduced. According to these investigations, the most important 
contribution to line blanketing stems from the iron group ele- 
ments. However, their abundances can be derived from ultravi- 
olet spectra only. The observed diversity of carbon and nitro- 
gen line strengths indicates that the treatment of metals has to 
be done carefully on a star-by-star basis. 

Apparently, the SPY sample suffers from observational se- 
lection biases. We found it likely that sdO stars were selected 
preferentially for the SPY target list. Hence sdO stars may 
be overrepresented with respect to sdB stars. To this end the 
Sloan Digital Sky Survey (SDSS) is a promising source of hot 
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subdwarf stars, as its selection criteria are very different from 
that of SPY. We ha ve already begun a spectroscopic analysis 
of SDSS sdO stars ( iHeber et all 12006 1) which will contribute 
significantly towards discriminating between the various evo- 
lutionary scenarios and hypotheses outlined above. 
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